Weak lensing mass distributions for 24 X-Ray Abell Clusters. 
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ABSTRACT 

We use the weak gravitational lensing effect to study the mass distribution and dynamical state 
of a sample of 24 X-ray luminous clusters of galaxies (0.05 < z < 0.31) observed with the FORS1 
instrument mounted on the VLT-Antu (Unit Telescope 1) under homogeneous sky conditions 
and subarsccond image quality. The galaxy shapes were measured in the combined V, I, R image 
after deconvolution with a locally determined point-spread-function, while the two-dimensional 
mass distributions of the clusters were computed using an algorithm based on the maximum 
entropy method. By comparing the mass and light distributions of the clusters in our sample, we 
find that their mass centers, for the majority of the clusters, is consistent with the positions of 
optical centers. We find that some clusters present significant mass substructures which generally 
have optical counterparts. At least in one cluster (Abell 1451), we detect a mass substructure 
without an obvious luminous counterpart. The radial profile of the shear of the clusters was 
fitted using circular and elliptical isothermal elliptical distributions, which allowed the finding 
of a strong correlation between the orientation of the major-axis of the matter distribution and 
the corresponding major-axes of the brightest cluster galaxy light-profiles. Estimates of how 
close to dynamical relaxation are these clusters were obtained through comparison of our weak- 
lensing mass measurements with the x-ray and velocity dispersion determinations available in the 
literature. We find that clusters with intra-cluster gas colder than 8 keV show a good agreement 
between the different mass determinations, but clusters with gas hotter than 8 keV present 
discrepant mass values. The clusters diagnosed to be out of equilibrium are Abell 1451, 2163 
and 2744, all of them having hints of substructure. Abell 2744 presents the largest discrepancy 
between its X-ray and weak-lensing temperature determinations, which can be interpreted as 
being due to the interaction between the two kinematical components along the line of sight 
found by Girardi & Mezzeti (2001). 

Subject headings: cosmology: observations- dark matter - galaxies: clusters: general- galaxies: clusters: 
individual (Abell 2744, Abell 1451, Abell 2163)- gravitational lensing 
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1. Introduction 

The emergence of mass structures in the uni- 
verse is currently linked to primordial density 
perturbations, dominated by cold dark matter, 
through a hierarchical clustering process involving 
gravitational instabilities. A natural consequence 
of this bottom-up scenario is that the most mas- 
sive structures form later in time, and, depend- 
ing on the world models, they are expected to 
show signatures of their assembling history. Since 
galaxy clusters are possibly the largest (nearly) 
virialized structures (e.g. Press & Schechter 1974), 
the study of their main properties (masses, shapes, 
radial profiles etc.) can provide valuable clues 
on the details of the agglomeration of matter in 
the universe and hence on the nature of the dark 
matter (e.g. Kauffmann et al. 1999). Clusters 
can also be used as cosmological probes provid- 
ing additional data that, together with cosmic mi- 
crowave background observations, can be used on 
the determination of several cosmological parame- 
ters. Studies of their mass function (e.g. Henry & 
Arnaud 1991; Rciprich & Bohringcr 2002), mass- 
to-light ratio (e.g. Bahcall Lubin,& Dorman; Carl- 
berg et al. 1996) and baryon fraction (e.g. White 
& Frenk 1991; Allen, Schmidt & Fabian 2002) can 
provide constraints on key cosmological parame- 
ters, like the mass density parameter, f^Mj the 
bias, and on the power-spectrum amplitude and 
shape parameter (Lacey & Cole 1994; Bahcall & 
Fan 1998). 

Galaxy clusters are complex systems that hold 
together galaxies, hot gas and dark matter. These 
components are governed by different physical 
mechanisms and their study requires the use of 
different observational techniques. Imaging in vis- 
ible light, reveals the cluster through their mem- 
ber galaxies which we know now that contribute 
only a small fraction of the cluster total mass. 
But, if these galaxies are in virial equilibrium, the 
depth of the cluster potential well can be accessed 
through their velocity dispersion, although the 
continuous accretion of field galaxies can bias the 
mass measurements to higher values (e.g. Sodre et 
al. 1989). The thin hot gas- mainly hydrogen- 
that permeates the cluster gravitational poten- 
tial is found to be at temperatures of the order 

1 Based on observations collected with ESO Very Large 
Telescope Antu (UT1) 



of 10 7 — 10 8 K, thus fully ionized and emitting 
X-rays via thermal bremsstrahlung process (see 
Sarazin 1988, for a more complete discussion). 
The emission in the X-ray band provides an effi- 
cient method to find galaxy clusters, and estimates 
of the cluster masses if the hypothesis that the gas 
is in hydrodynamical equilibrium is adopted. The 
advantage of using the hot intracluster gas rather 
than the member galaxies for mass determinations 
is that it has a much shorter relaxation time due 
to its self-interactivity (Sarazin 1988). Neverthe- 
less, there is evidence, that many clusters present 
significant departures from dynamical relaxation 
(e.g. Geller & Beers 1982; Sodre et al. 2001). Pre- 
sumably the clusters that are in the top or the 
mass function will achieve last their formation pro- 
cess, so we expect to find some degree of correla- 
tion between the departures from equilibrium and 
the mass of the clusters. 

The development of gravitational lensing tech- 
niques, for both the strong and weak regimes, has 
presented a new way of measuring masses regard- 
less the nature or dynamical state of the matter 
(Fort & Mellier 1994; Mellier 1999). Strong lens- 
ing, that relies on the modeling of systems of mul- 
tiply lensed background galaxies, has proved to be 
an accurate method, mainly when the number of 
gravitational arcs is large (e.g. Kneib et al. 1996). 
However, this technique can probe only the very 
central regions of the clusters, and is limited to the 
most massive and concentrated objects. On the 
other hand, the statistical analysis of the weak dis- 
tortion caused by the cluster shear field on images 
of faint background galaxies allows the mapping of 
the matter distribution to much larger radius; be- 
sides, this gravitational distortion is detectable in 
almost all clusters when large telescopes with cur- 
rent instrumentation are used. After the seminal 
work of Tyson, Valdes & Wenk (1990), the weak- 
lensing technique started to be used for the study 
of some clusters having special interest because of 
extreme values in some of their measured proper- 
ties, and new methodologies were established (see 
Mellier 1999, for a review). 

With the development of wide field cameras and 
the increase in the number of large telescopes, it is 
now possible to study large samples of galaxy clus- 
ters selected using well defined criteria, ideally ob- 
served under similar conditions, and analyzed with 
the same technique. Dahle et al. (2002) studied 38 
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Northern X-ray luminous clusters using data col- 
lected with 2-meter class telescopes and selected to 
be representative of the most massive clusters in 
the 0.15 < z < 0.35 redshift range. However, the 
clusters in this sample were not all observed with 
the same instruments and, sometime, with varying 
atmospheric conditions. In this paper we present a 
weak-lensing analysis of a sample of 24 X-ray lu- 
minous southern clusters with 0.05 < z < 0.31, 
based on imaging observations with the ESO- 
VLT telescope taken under sub-arcsecond condi- 
tions. The clusters were selected from the X-ray 
Bright Abell Cluster Survey (XBACS Ebeling et 
al. 1996), with -50° < 6 < 15°, 12 h < a < 1 h, 
Lx > 5 x 10 44 erg s . Similar observations have 
been done of a complementary sample of 27 X-ray 
clusters, with 1 h < a < 12 h , whose weak-lensing 
analysis is underway. Hereafter, and due to the 
rather rare choice of a lower-z limit of 0.05, we will 
refer to our whole VLT survey as the Low-z Lens- 
ing Survey of X-ray Luminous Clusters (LZLS), 
and the sample presented here is its Part I. The 
data used in this paper derive from observations 
primarily designed for detection of strong-lensing 
features but conceived for weak-lensing measure- 
ments also; the service observing mode at ESO was 
essential to insure the high quality and homogene- 
ity of the images. The statistics of the occurrence 
of bright arcs for the whole LZLS survey will be 
presented in a separate paper (Campusano et al. 
2003), together with its implications with respect 
to the clusters inner mass radial profile. In future 
papers, wc will address the combined evaluation of 
the mass distribution we will use both weak and 
strong lensing in those LZLS clusters with gravi- 
tational arcs focusing on those involving arcs with 
redshift determinations. 

In this paper we determine the mass distribu- 
tions for the galaxy clusters in our sample and 
their total masses using weak-lensing techniques, 
and investigate their dynamical state through 
comparison of the weak-lensing masses with al- 
ready published virial and X-ray mass estimates. 
In Section 2 we describe the sample selection and 
the observations. In Section 3 we present the 
procedures adopted in the weak-lensing analysis, 
including galaxy shape measurements and the re- 
construction of cluster density maps, as well as 
the results obtained. The discussion of the results 
is presented in Section 4 together with the com- 



parison with dynamical, X-ray and weak-lensing 
masses taken from the literature. In Section 5, 
we summarize our main conclusions. In the Ap- 
pendix, we display images, mass and light maps, 
and weak-shear profiles for the cluster sample. 

Throughout this paper we adopt Q,m = 0.3, J7a 
= 0.7 and H = 70 km s" 1 Mpc" 1 . 

2. Observations 

2.1. Sample Selection 

The galaxy clusters investigated in this work 
were selected based on their high X-ray luminosi- 
ties (Lx). Initially, wc selected clusters in the 
X-Ray Brightest Abell Cluster Catalog (XBAC) 
(Ebeling et al. 1996) with L x > 5 x 10 44 erg s" 1 
in the 0.1-2.4 keV band, which constitutes a good 
threshold to identify real, massive clusters (e.g. 
Luppino et al. 1999). We then restricted the sam- 
ple to clusters that could be reached with ESO 
VLT (-50° < S < +15°). A low-z limit of 0.05 
was chosen, because for lower redshifts the angular 
size of the clusters is too large compared with the 
field of the imaging camera. For z < 0.05, possible 
gravitational arcs would appear projected deep in- 
side the central galaxy image and the angular size 
of cluster galaxies in the camera field would be 
very high, making hard to find a suitable num- 
ber of background galaxies (without strong light 
contamination) for a weak-lensing analysis. ESO 
VLT was considered a very good choice because 
the observations could be completed with some 
10 hours of service observing and the homogene- 
ity and image quality of the observations could be 
optimally achieved. These criteria produced 27 
targets within 12 h < a < 1 li. Three of these 
clusters (A68, A1689 and A1835) have been al- 
ready observed with HST and CFHT by one of us 
(JPK) and the results will be published elsewhere 
(Smith et al. 2003; Bardeau et al. 2003). The re- 
maining 24 rich galaxy clusters that are studied in 
this paper are presented in Table 1. This sample 
constitutes the Part I of the VLT Low-z Lensing 
Survey of X-ray Luminous Clusters (LZLS). 

2.2. Imaging and Data Reduction 

All clusters were observed in service mode us- 
ing VLT Antu (Unit Telescope 1), from April to 
July 2001 (ESO program 67.A-0597), with the 
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Table 1 

VLT Massive Abell Cluster Sample 
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Note. — (1) Cluster name. (2), (3) Equatorial J2000 coordinates. (4) Rcdshift. (5) X-ray 
luminosities in the [0.1-2.4] keV band (Ebeling et al. 1996). (6) FWHM of stellar images in 
the combined VRI image. (7) Surface brightness corresponding to la above the sky average in 
the R image. (*) A2744 is also known as AC118. 
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F0RS1 instrument (FOcal Reducer/low disper- 
sion Spectrograph) working in its imaging mode. 
FORS1 employs a TK2048EB4-1 backside illumi- 
nated thinned CCD, with 2048 x 2048 pixels, each 
with 24 x 24/x 2 area. The employed imaging mode 
of FORS1 results in a pixel scale of 0".2 and a 
square field of view of 6'. 8 on a side, correspond- 
ing to 0.4, 0.9 and 1.8 Mpc at the smaller, me- 
dian and greater redshifts of the clusters of our 
sample, which are 0.05, 0.126 and 0.305 respec- 
tively. The imaging was chosen to be centered on 
the clusters cores and done through the V, R, and 
I-bands with exposure times of 330s in each filter, 
ensuring a good detection of faint galaxies in all 
bands. The color information is essential to dis- 
criminate elongated objects belonging to the clus- 
ter from the background sources needed for the 
weak-lensing measurements. All the observations 
were conducted under good sky transparency and 
and an excellent image quality characterized by a 
median seeing of 0".6, ranging from 0".4 to 1".0. 

The data reduction were done in a standard way 
using IRAF packages. For each cluster we com- 
bined the images in the three filters, normalizing 
by their modes, to produce a high signal-to- noise 
"VRI" image that was used for the weak-lensing 
analysis. The images in each band were used for 
the photometry of the objects in each cluster field. 
Comparing stellar shapes measured in the VRI im- 
ages with the single-band ones, we concluded that 
the combination process did not introduce biases. 
Tests on galaxies revealed that the use of the com- 
bined VRI image improves the precision in the 
shape measurements in comparison to single-band 
measurements. 

Zero-point magnitudes, extinction, and color 
coefficients for transformations from instrumental 
to standard magnitudes were provided by the ESO 
pipeline. Almost all data have been collected un- 
der photometric sky conditions or nearly so. Ta- 
ble 1 contains, besides general informations about 
each cluster in our sample, the seeing (FWHM) 
derived from stellar images in the combined VRI 
image. The table also includes the surface bright- 
ness corresponding to a ler fluctuation above the 
sky background in the R image, which is a measure 
of the quality and depth of these images. 



3. Weak-Lensing Analysis 

3.1. Construction of the Catalog 

The detection of objects in the images and the 
extraction of their main parameters were made 
with SExtractor (Bertin & Arnouts 1996). We 
run this program separately for the single band 
images and for the combined one. The four cata- 
logs obtained were then combined by matching the 
Cartesian position of the objects, with a tolerance 
of 0".6. 

Magnitudes (Sextractor's MAG_BEST) were 
measured in the single band images, and the colors 
were computed using 4" diameter apertures. As- 
trometric and morphological data were obtained 
from the VRI image. 

We have produced three sub-samples of de- 
tected objects: stars, cluster galaxies, and other 
galaxies. Stars have been selected using two dif- 
ferent criteria. For unsaturated objects (R^ 20.0 
mag) we have used the FWHM of their light pro- 
files. In Figure 1 we show an example of the distri- 
bution of this parameter as a function of R mag- 
nitude, which shows that stars can be easily sepa- 
rated from galaxies down to very faint magnitudes. 
For brighter objects, we have used the SExtrac- 
tor neural network classifier. 

In the absence of spectroscopic redshifts for 
most of the galaxies in the clusters, the cluster 
galaxies were selected from the color-magnitude 
diagram (V — 7) X I for each cluster (Figure 2). In 
this type of diagram, the early-type cluster galax- 
ies occupy a well defined locus, the so called red 
cluster sequence, with almost the same color. Af- 
ter determining the typical color of the ellipti- 
cals in each cluster, we selected as cluster galax- 
ies those with colors within a strip of width 0.2 
mag matching the red sequence of each cluster. 
This selection procedure naturally will not include 
the bluer galaxies belonging to a particular clus- 
ter and having colors outside the color range de- 
fined by the strip in the cluster color-magnitude 
diagram. These bluer galaxies not recognized as 
cluster members cannot be distinguished from the 
background galaxies needed for the weak-lensing 
analysis and thus may introduce a certain level of 
"noise" . We will return to this issue later. 
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Fig. 1. — Magnitude - FWHM diagram for objects 
in the field of A 1451. Filled squares represent ob- 
jects classified as point-like sources, following the 
criteria described in the text. The continuous line 
is the threshold FWHM adopted for star/galaxy 
separation of objects fainter than R=20. The star 
sequence for unsaturated objects is clearly identi- 
fiable below the line. 

3.2. Shape Measurements 

Wcak-lensing analysis needs careful measure- 
ments of the ellipticity of background galaxies ob- 
served in the field of a cluster of galaxies. 

We made galaxy shape measurements, includ- 
ing the correction of the seeing circularisation and 
PSF anisotropies, using the im2shape software 
(Bridle et al. 2002). Following a proposition by 
Kuijken (1999), this software models astronomi- 
cal objects with a sum of Gaussians with elliptical 
base. Let x be a position on a CCD image. The 
intensity of a galaxy as a function of position is 
modeled as 

/ w = E^^ (x_Xi)TC!(x_Xi)/2 « 

. 27T Gj 

2 ' 1 

where the Xj, A4, and Cj are, respectively, the cen- 
ter of each Gaussian, its amplitude, and its covari- 
ance matrix. The components of C, can be writ- 
ten is terms of the ellipse parameters a, b, and 8, 
that are its semi-major and semi-minor axes and 
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Fig. 2. — Color-magnitude diagram for galaxies in 
the field of A1451. Galaxies belonging to the red 
cluster sequence are shown as filled squares. The 
vertical lines shows the range of magnitude for the 
galaxies used for the lensing analysis. 

the position angle of the major axis, respectively. 
The ellipticity, following the convention usually 
adopted in lensing studies, is defined as 



At this point it is convenient to define the el- 
lipticity projected tangentially to the direction of 
a predefined cluster mass center, 

et = e cos [2(0-0)] (3) 

where <f> is the azimuthal coordinate of the galaxy. 

im2shape uses a Bayesian approach to deter- 
mine these parameters, with carefully calculated 
uncertainty estimates. It also deconvolves the 
measured shape using a PSF that is also given 
as a sum of Gaussians, so that the dcconvolution 
process can then be done in a fully analytical way. 

The steps actually performed in measuring 
galaxy shapes on our images can be summarized as 
follows. First, im2shape is used to model images 
of the unsaturated stellar objects as single Gaus- 
sians. These stars (in fact stellar-like objects) 
are divided in four sets, following their positions 
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regarding the CCD quadrants. Stars with e, 
or base area unrepresentative of their respective 
quadrant averages are removed through a 2-sigma 
clipping process. This cleaned star catalog is used 
to map the PSF along the field. After, the same 
program is applied to the galaxies in the CCD im- 
age. Each galaxy image is modeled as a sum of two 
Gaussians, which is adequate for exponential and 
de Vaucouleurs profiles, given typical galaxy sizes, 
CCD pixel-sizes and PSFs (see Kuijken 1999, for 
a discussion). The parameters x, y, e, and 9 are 
forced to be the same for both Gaussians. Using 
the five nearest stars, a local PSF is calculated for 
the position of each galaxy and a deconvolution 
of their images is performed. The result of this 
procedure is a list with the elliptical parameters 
of the galaxies in a cluster field. 

3.3. Weak-Lensing Sample Selection 

A key step in weak-lcnsing analysis is the se- 
lection of galaxies that will be used as probes of 
the gravitational field of the lensing cluster. These 
galaxies must be behind the cluster and we should 
be able to measure their shapes with good accu- 
racy. For this, it is necessary to have distance 
estimates for all objects in the field. It is possible 
to select background galaxies through photometric 
redshifts (e.g. Athreya et al. 2002), but an efficient 
application of this technique requires more abun- 
dant color information than what we have. In the 
absence of redshift information, we proceeded by 
building a master catalog of putative background 
galaxies, selected by flux in the non-cluster galaxy 
catalog. We selected only galaxies with absolute 
R-band magnitude fainter than -16.6 at the cluster 
distance (roughly M*+5; Goto et al 2002). 

The final catalog is obtained after removal of 
a number of objects, obeying the following crite- 
ria. Initially, all objects inside a radius of 1 arcmin 
around the cluster center were excluded from the 
lensing sample (following, for example, Hoekstra 
et al. 2002). The reason for this is twofold. First, 
in these regions (mainly the inner 30") mass den- 
sities are closer to the critical value, and then the 
wcak-lensing regime is no longer a good approxi- 
mation. Second, in central regions there is a high 
surface density of low luminosity cluster compo- 
nents and diffuse light, that can contaminate the 
wcak-lensing sample and bias the shape measure- 
ments. 



Relevant cases of light contamination by lumi- 
nous neighbors were detected by comparing dif- 
ferences in the object center estimated by S Ex- 
tractor and by IM2SHAPE, because the first has 
a deblending procedure, while the second does not. 
Thus, objects with differences between center es- 
timates greater than 3 pixels were also removed 
from the sample. Finally, all galaxies with uncer- 
tainties in e t larger than 0.35 were also removed. 
We found that this value is a good compromise 
between the total number of galaxies in the final 
sample and the accuracy of shape measurements. 
It should be stressed, however, that the results 
presented here depends only weakly on this value, 
since in the whole analysis the inverse of the square 
of the uncertainties are used as statistical weights 
(see for example Section 3.5). 

The total number of galaxies useful for shear 
measurements is presented in column (2) of Ta- 
ble 2. Their surface density ranges from 4.7 to 
17.8 gal/arcmin 2 . This broad range are due to 
cosmic variance, almost unavoidable differences in 
the deepness of the images, (see the values of the 
seeing and fj- s k y in Table 1) and also to the dif- 
ferent bright apparent magnitude cutoffs for the 
galaxy lensing samples. 

3.4. Mass and Light Maps 

In the weak-lcnsing regime, the mass distribu- 
tion can be obtained from the pattern of distortion 
of galaxy shapes at a position x, (e t (x)), that de- 
pends only of the reduced shear g = — k), 
where 7 and n are the shear modulus and conver- 
gence at x. The latter is proportional to the mass 
density and is defined as: 

where the critical surface density S c is given by 
c 2 D, 

and Di, D s , and D; s are the angular diameter 
distances from the observer to the lens, from the 
observer to the source, and from the lens to the 
source. 

We have reconstructed the two-dimensional 
mass density distributions from distortion maps 
using the LensEnt software (Bridle et al. 1998; 
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Marshall et al. 2002). This software computes k 
maps using a maximum-entropy method, taking 
each background galaxy image shape as an in- 
dependent estimator of the reduced shear field. 
This map is smoothed by a Gaussian function, 
called intrinsic correlation function (ICF), whose 
FWHM is determined by Bayesian methods. The 
program avoids any binning of the data and allows 
the reconstruction of complex mass distributions. 
We present in column (3) of Table 2 the optimal 
ICF-FWHM for the smoothing function found by 
LcnsEnt. 

The resultant mass maps for all clusters, with 
64 x 64 pixels across, can be seen in Appendix A. 
In order to minimize edge effects, LcnsEnt pro- 
duces maps in a grid twice as big as the data re- 
gion. Here we only show the part of a map corre- 
sponding to the field actually observed. 

It should be noted that since this method takes 
into account only galaxy distortions and ignores 
their amplifications, the convergence determined 
this way will suffer from the "mass sheet degener- 
acy" : the distortion pattern can be reproduced by 
any mass distributions that obeys the transforma- 
tion k — > Xk + (1 — A), where A is a real number. 

To estimate the significance of the features 
found in mass maps we have also performed 100 
bootstrap realizations for each cluster, adding to 
galaxy ellipticities Gaussian noise with amplitude 
equal to their observational errors. From the mass 
maps produced in the simulations, we computed 
a standard deviation (er) map and, by dividing 
the mass map by the <r-map, we obtained a sig- 
nificance map. The significance maps show that 
most clusters presents individual pixel signal-to- 
noise ratios ranging between 1 and 3, with values 
as high as 8 being found in some cases. 

We have also computed light maps (also shown 
in Appendix A) of the light distribution of galaxies 
in the red sequence. The field is divided in 32 x 32 
pixels and the smoothed light density of each pixel 
is estimated as 



C cx Li 



(6) 



where Li is the luminosity of the i-th galaxy, and 
di is its angular distance to the pixel center. The 
value of a is chosen to have the same FWHM of 
the ICF of the mass map. 



For two clusters, A1651 and A1664, the algo- 
rithm was not able to reconstruct the density dis- 
tribution. This will be discussed in Section 4.5. 




Fig. 3. — Mass (left) and Light (right) maps for 
the cluster A1451. The mass map is produced us- 
ing the distortions of the faint background galaxies 
as probes of the mass density of the cluster. The 
gray level scale shown at the right edge of the map 
is in unities of the critical density. The light map 
is the result of the smoothing of the magnitude 
of the clusters members belonging to the red se- 
quence (see text for details). 



3.5. Profile Fitting 

In order to avoid the mass sheet degeneracy and 
to obtain quantities that can be compared with 
data determined by other methods, wc fit physi- 
cally motivated parametric models to the distor- 
tion maps. We considered tho different mass mod- 
els: a singular isothermal sphere (SIS) and a sin- 
gular isothermal ellipsoid (SIE). 

Let us assume a polar reference system with 
origin at the cluster center, so that any point on 
the image can be represented by an angular radial 
coordinate 6 and an azimuthal angle <j>. For the 
SIS profile the convergence and the shear are given 

by 

where 9e is the Einstein radius, which is related 
to the one-dimensional velocity dispersion of the 
isothermal sphere (a sis) as 



4tt 



j sis 
c 2 D s 



(8) 



The SIE profile (e.g., Kormann, Schneider & 
Bartelmann 1994) has properties similar to those 
of SIS models, but the values of the shear and the 
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convergence are now given by 

K = 7 = TW f [cos ^ ~ a) + f sin {<t> ~ a)r V2 

(9) 

where / < 1 is the axial ratio b/a, a is the po- 
sition angle of the ellipsoidal matter distribution 
and <f) the azimuthal coordinate. For this model, 
in the weak-lcnsing regime (re << 1), the shear is 
oriented tangcntially to the direction to the mass 
center. 

The SIS model has only one free parameter, Be, 
whereas the SIE has three free parameters: 9e, /, 
and a. In both cases the position of the cluster 
center were defined in advance with the help of 
the mass maps. For each model, the best-fit pa- 
rameters were obtained through minimization of 
the x 2 statistic, defined as: 

where et,i and <7 et ,i are the tangential ellipticity 
and its error for the i-th galaxy, <r e = 0.3 (Mel- 
lier 1999) is the dispersion associated with galaxy 
intrinsic shapes, and g t ,i is the tangential reduced 
shear at the position of this galaxy, which quan- 
tifies the ellipticity induced by the lensing distor- 
tion. For all models we consider here it may be 
assumed that, in the weak lensing regime, \g\ = g t 
This method of profile fitting is similar to the one 
developed by King & Schneider (2001) and has the 
advantage of avoiding data binning. 

The mass center was defined in general as the 
center of the mass map. When there is consistency 
between the peak of the mass map and a bright 
cluster galaxy (the dominant one in most cases), 
the latter one was used as the center position. In 
some cases (A2104, A3739 and A3888) the maxi- 
mum of the mass map and therefore the adopted 
mass center, is in a position between bright clus- 
ter members, and in other cases (A85, A2811 and 
A1650) in a position uncorrelated with the cluster 
galaxies. 

The best-fit parameters of SIS and SIE mod- 
els are also shown in Table 2. We have estimated 
&SIS {vsie for the elliptical model) from 9e us- 
ing for each cluster an average value of Di s /D s 
obtained in a way analogous to that adopted by 
Dahlc et al. (2002). Considering the minimum 
and average values of I magnitudes in our weak- 
lensing sample, we selected sub-samples of HDF 



galaxies (Fcrnandez-Soto, Lanzetta & Yahil 1999) 
and, from their photometric rcdshifts, the mean 
value of the distance ratio above was calculated. 
In this table we also show the average redshift 
of the galaxies used in the lensing analysis with 
their respective uncertainties (that are not taken 
into account on the estimate of the uncertainties 
in (t sis ). 

When fitting noisy data or clump mass distri- 
butions with the SIE profile, one often obtain un- 
realistic low values of the axial ratio (/ < 0.3). 
Clusters where this happened have an empty en- 
try in Table 2. We dealt with it adopting a con- 
stant value of / = 0.3 when the fit gives values 
smaller than this one. On the other hand, the 
value obtained for the position angle is quite ro- 
bust, being rather insensitive to /. Soft core (i.e. 
non singular) and Navarro, Frenk & White (1997) 
profiles have parameters, such as the core radius or 
the NFW concentration parameter, whose deter- 
mination depends critically on information from 
the inner parts of the cluster. Since we do not 
have useful data in the central arcmin, or even at 
larger radius, (that are important for NFW pro- 
files) , we are unable to fit properly our data using 
these models. 

Another point that deserves mention in Table 
2 are the small values of xf-ed- This happens due 
to the large value of the intrinsic shape dispersion 
cr e . 

4. Discussion 

4.1. Mass and Light Distributions 

Although each galaxy cluster of our sample has 
particular features, as can be seen in the display 
of the mass and light distributions provided in 
the Appendix, there are some general trends and 
correlations that can be drawn from these maps. 
A detailed and more quantitative comparison of 
these distributions will be addressed in a future 
paper. 

In most cases (17 out of 22), the center of mass 
of the cluster (nearly) coincides with the position 
of the brightest cluster galaxy (BCG) or the main 
galaxy clump. We are including in this category 
the clusters A22, A1437, and A2163, where the 
mass maps show peaks near the position of the 
central galaxy, despite the presence of more promi- 
nent peaks near the edges of the frame (see rele- 
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Table 2 
Best fit parameters 











SIS 






SIE 






Cluster 


N ICF-PWHM 


<*> 


Oe 


°sis 




O'SIE 


/ 


Q 


Xred 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



A2744 


323 


220 


1 


,01±0 


,04 


27.5±4.3 


1491 


± 


116 


1 


.54 


1545 


± 


111 


0.52 ±0.21 


1 


.73 


± 


16 


.65 


1 


.53 


A22 


554 


240 





,84±0 


.05 


10.7±3.6 


663 


± 


113 


1 


.46 


812 


± 


79 




14 


81 


± 


10 


.79 


1 


.43 


A85 


492 


110 





81±0 


.05 


20.7±3.8 


917 


± 


85 


1 


11 


926 


± 


82 


0.54 ±0.27 


21 


37 


± 


21 


80 


1 


11 


A2811 


498 


80 





,86±0 


.05 


16.3±5.6 


863 


± 


149 


1 


.13 


861 


± 


134 




-34, 


38 


± 


15 


22 


1 


13 


A1437 


385 


160 





,84±0 


,05 


10.8±4.5 


734 


± 


152 





.95 


780 


± 


127 




45, 


94 


± 


16 


74 





.94 


A1451 


462 


80 





,90±0 


,05 


16.8±3.9 


993 


± 


114 





.95 


1087 


± 


98 




-14, 


60 


± 


9 


80 





.94 


A1553 


178 


160 





81±0 


,05 


15.4±6.2 


923 


± 


186 


1 


.05 


925 


± 


164 




-49, 


04 


± 


19 


70 


1 


.05 


A1650 


279 


160 





,80±0 


,05 


17.7±5.4 


876 


± 


133 


1 


.08 


1033 


± 


109 




-29, 


03 


± 


10 


.36 


1 


.06 


A1651 


391 







,73±0 


,04 


1.0±4.9 


210 


± 


498 





.91 


489 


± 


116 




70 


72 


± 


47 


.99 





.91 


A1664 


246 







,80±0 


,05 


-1.9±6.3 












402 


± 


264 




-74, 


14 


± 


92 


74 


1 


.05 


A2029 


549 


140 





,83±0 


,05 


25.2±3.1 


1039 


± 


64 





.86 


1052 


± 


63 


0.62 ±0.21 


-26, 


11 


± 


19 


.63 





.86 


A2104 


378 


90 





,86±0 


,05 


20.7±4.3 


1039 


± 


108 





,86 


1038 


± 


102 


0.39 ±0.22 


58, 


44 


± 


14 


44 





.86 


A2163 


261 


240 





,89±0 


,04 


17.6±5.0 


1021 


± 


146 





.90 


1094 


± 


125 




-83, 


53 


± 


12 


50 





.88 


A2204 


347 


110 





,86±0 


,05 


20.3±4.1 


1028 


± 


104 





.81 


1035 


± 


102 


0.58 ±0.32 


9 


.73 


± 


27 


80 





,81 


A3695 


294 


130 





,83±0 


,05 


19.5±4.9 


928 


± 


117 


1 


.10 


987 


± 


104 




-23, 


08 


± 


11 


08 


1 


.09 


A3739 


229 


110 





,89±0 


,05 


26.2±4.9 


1203 


± 


112 


1 


.06 


1210 


± 


96 




31, 


71 


± 


7 


.37 


1 


.05 


A2345 


364 


150 





,89±0 


,05 


15.0±4.6 


909 


± 


138 


1 


.08 


965 


± 


126 




87, 


18 


± 


14 


14 


1 


.08 


A2384 


420 


150 





,80±0 


,05 


12.2±4.2 


737 


± 


126 


1 


.05 


787 


± 


108 




18 


97 


± 


14 


18 


1 


,05 


A2426 


343 


210 





,80±0 


,05 


10.1±4.8 


676 


± 


158 





.91 


767 


± 


132 




58, 


74 


± 


18 


14 





.90 


A3856 


429 


110 





,87±0 


,05 


20.1±4.0 


1005 


± 


97 


1 


.04 


1022 


± 


93 




-33 


40 


± 


8, 


,86 


1 


.03 


A3888 


270 


130 





,86±0 


,05 


18.6±4.7 


981 


± 


125 


1 


.21 


1008 


± 


125 


0.58 ± 0.40 


79, 


13 


± 


34 


,82 


1 


.22 


A3984 


607 


210 





,95±0 


,05 


22.2±3.1 


1093 


± 


77 


1 


.27 


1150 


± 


66 




55. 


19 


± 


6 


.38 


1 


.25 


A2597 


524 


150 





,89±0 


,04 


14.1±3.7 


776 


± 


101 


1 


.03 


853 


± 


79 




-18, 


57 


± 


10 


.59 


1 


.01 


A4010 


673 


160 





,87±0 


,04 


11.3±3.3 


706 


± 


104 


1 


27 


776 


± 


84 




56, 


08 


± 


11 


.95 


1 


.26 



Note. — (1) Cluster name. (2) Number of galaxies used for weak-lensing analysis. (3) FWHM of the smoothing 
function used to build the mass map in arcsec. (4) Average redshift and standard deviation of the galaxies in the weak 
lensing sample (see text for details). (5,6) The Einstein radius (in arcsec) and asis and respective uncertainty obtained 
by the fitting of a SIS profile. (7) The \ 2 of the SIS fit divided by the number of degrees of freedom. (8), (9), (10) 
<?sie, axial ratio and position angle (counted counterclockwise with the North direction been the origin) and respective 
uncertainties obtained by the fitting of a SIE profile. When the best SIE fit gives a value of f smaller than 0.3, we have 
done the fit again using a fixed value of f=0.3 and let the entry for f in the table empty. (11) The \ 2 of the SIE fit 
divided by the number of degrees of freedom. 
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B (h 70 1 kpc) 
400 600 



C7 SIS = 993±114 km s 




6 (arcmin) 

Fig. 4. — The shear profile for the cluster A1451. 
The filled symbols correspond to average ellipticics 
of the faint background galaxies projected tangen- 
tially to the cluster center. Each point represents 
nearly 1/5 of the galaxies. The open squares are 
the same, but of ellipticies projected in a direction 
45°relative to the center. The solid line shows the 
best fit SIS profile. 

vant figures in the Appendix). 

Although most of the clusters mass maps show 
single clumps, an effect that can be attributed in 
part to poor resolution, some of them present sec- 
ondary peaks (e.g. A1451, A3695), or extensions 
of the main peak in the direction of substructures 
present in the light distribution (A2744, A2104, 
A2597 and A3739). Conversely, the clusters A1553 
and A2163 present prominent light clumps with- 
out counterparts in the mass map. 

For some cluster, the mass and light maps sug- 
gest the existence of "dark clumps" (mass concen- 
trations without optical counterparts), but their 
reality is doubtful since most of them are near the 
edge of the observed cluster field. Nevertheless, in 
the map of A1451 (see Fig. 13) there is a signif- 
icant structure, South of the main cluster clump, 
with an extension to the West, that is not too close 
to the frame border and is not associated with 
bright cluster galaxies. Indeed, our bootstrapping 
of galaxy ellipticities shows that this is a robust 
feature. The X-ray isophotes of Valtchanov et al. 
(2002) also do not show any emission excess at 



this region. The nature of this feature, whether 
it is a real dark clump belonging to A1451 or a 
background massive structure, will be subject of 
future investigations. 

Dark matter profiles tend to present significant 
departures from circular symmetry. This can be 
noticed by the smaller values of Xred presented by 
SIE in relation to the SIS fits, as well as by the 
uncertainties of ctsie, that are ^20% smaller for 
the non-circular profile. Furthermore, the mass 
maps often show elongated distributions, with ma- 
jor axis with position angles similar to those given 
by the SIE fitting. Sometimes the ellipticity in- 
ferred from the mass map is smaller than that from 
the SIE fitted profile. This is due to two reasons: 
first, the smoothing process in the mass map cal- 
culation tends to circularize features whose scales 
are smaller than the map resolution; second, as 
discussed before there is a bias toward small el- 
lipticity values in the fitting of SIE profiles with 
noisy data. 

A very robust result that comes from the SIE 
fits is that the major-axis of the mass distribution 
is clearly aligned with the major-axis of BCGs. 
This is shown in Fig. 5, where we present the 
absolute value of the difference between the posi- 
tion angle of a BCG (from SExtractor)) and the 
position angle given by the SIE mass model. For 
clusters without a clear dominant galaxy (A2744, 
A1451, A3888, A3984), the direction linking the 
two first ranked cluster members was adopted to 
characterize the orientation of the light distribu- 
tion. On the other hand, for at least two cases 
(A1451, A2426) the main direction of the dark 
matter distribution is more related with the di- 
rection linking bright cluster members than with 
the direction of the BCG. 

The strong alignment between light and mass 
in Mpc scales found here has also been detected, 
at smaller scales, by the excess of strong lensing 
features found in the direction of the major axis 
of BCGs (Luppino et al. 1999; Campusano et al. 
2003). 

Considering that alignments between several 
cluster components (BCGs - cluster galaxies - 
sub-structures - surrounding large-scale struc- 
tures) have already been found (e.g. West, Jones 
& Forman 1995; Plionis & Basilakos 2002), the 
present result reinforces the conclusion that clus- 
ters, in the same way as their dominant galaxies, 
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Fig. 5. — Alignment between the position angle 
of the brightest cluster galaxy asca and of the 
mass distribution asiE- For clusters without a 
clear dominant galaxy, the direction linking the 
two first ranked cluster members is used instead. 

have been formed by the anisotropic merger of 
smaller structures along large-scale filaments, in 
close agreement with numerous cosmological sim- 
ulations based on the CDM paradigm (e.g. Jenkins 
ct al. 1998). It should be stressed here that the 
diffuse luminous halo associated with a cD galaxy 
tends to produce the opposite effect. This halo will 
produce a significant gradient in the background 
of a certain galaxy, that contaminates the measure 
of its shape, causing a stretch of it in the radial 
direction. As this halo is brighter in the direction 
of the major-axis, galaxies located in this direc- 
tion will be more affected. This is probably what 
is happening in A2029, which has a very luminous 
diffuse component (Uson, Bough & Kuhn 1991) 
and where the major-axis of the mass distribution 
is about 60 degrees of the central cD major-axis. 
The hypothesis of light contamination in this case 
is supported by the fact that when the central 
exclusion radius is increased from 60" to 75", 
it produces a major change in a, that becomes 
almost parallel to the BGC main axis, however 
without causing appreciable changes in ctsie- We 
did not detect such effect in any other cluster of 
the sample. 



4.2. The dynamical state of massive clus- 
ters 

The dynamical state of the clusters in the 
present sample can be accessed through a compar- 
ison between the mass estimates made by wcak- 
lensing, dynamical, and X-ray techniques. The 
key idea is that the weak-lensing signal is inde- 
pendent of the cluster dynamical state, whereas 
other techniques rely on assumptions of equilib- 
rium. 

We have collected measurements of Tx and a v 
from the literature (see Table 3) . Assuming energy 
cquipartition between cluster galaxies and gas, we 
may obtain an estimate of the temperature from 
our weak-lensing determination of asis (and for 
o'Sie), through the relation (e.g. Sarazin 1988) 



where fi = 0.61 is the mean molecular weight and 
ma the hydrogen mass. The results are also in Ta- 
ble 3. In Fig. 6 we present a comparison between 
temperatures computed in this way with those ac- 
tually obtained with X-ray observations. In Fig. 
7 we compare the velocity dispersions determined 
through lensing with those obtained directly from 
galaxy velocities. The first point to note in these 
figures is that for many clusters the weak-lensing 
results agree with X-ray or dynamical data within 
1.5 ct. The clusters A1451, A1553, and A2163, 
however, show temperatures and velocity disper- 
sions significantly in excess when compared with 
estimates derived from weak-lensing data, sug- 
gesting that they should be dynamically active. 
Actually, the analysis of the A2163 temperature 
map made by Markevitch & Vikhlinin (2001) with 
Chandra data, shows at least two shocked regions 
and other evidences that the central region of this 
cluster is in a state of violent motion. In the same 
way, Valtchanov et al. (2002) describes A1451 as 
being in the final stage of establishing equilibrium 
after a merger event, whereas its high X-ray tem- 
perature (13.4 keV) would be probably due to a 
shock occurred recently. A2744 seems to be an ex- 
ception, since its temperature is significantly lower 
than the weak-lensing estimate. The nature of this 
cluster will be discussed in more detail in Section 
4.4. 

We present in Table 4 the weighted (by the 
error of the ratio) mean ratios Tx /Tsis,sie and 
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Table 3 



Nome 


kTsis 
(keV) 


(keV) 


(keV) 


Refs. 


a v 
(km s- 1 ) 


N 


Refs. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


( 7 ) 


(8) 


A2744 


1 a 1 o -I - 2 . 28 

14.13I 211 


1 r 1 t+2. 26 

15.17I 210 


i i n/i+0.49 


(a) 


1777 -125 


■55 


(g) 


A OO 
AZZ 


9 7 q + 1.03 
z - ' y -0.87 


/I Ofi+0.78 
4 -26_o.73 






t?no i qc 1 

D9o±zol 


1 


fv, \ 

w 


A85 


r O/I+1.04 

°-' 34 -0.94 


r A r+1.00 

°- 40 -0.92 


a r i +0.16 
D - O1 -0.23 


(b) 


1097I 63 


305 


(1) 


A OO 1 1 

A2811 


4.73I 149 


4 - 76 -1.33 


r qi +0.17 
°- ,il -0.16 


(c) 


cnc+200 


1 Q 

13 


u) 


A 1 A 0*7 

A143/ 


3.42I 127 


o O7+1-30 

3.87_ 1 n 












A1451 


6 - 27 -1.36 


r, ^+1.40 

7 - 47 -1.28 


13.40]^ 5 


(a) 




3 1 


(a) 


A1553 


r ,io+2.40 

5-421!. 96 


r ■) — h2.00 
-1.69 


n i f + 1.02 
9-16_o.64 


(c) 








A1650 


4.<3<5_ 1 37 


7 r+1.54 

' - zo -1.39 


r RQ + 0.30 

5.68_ 027 










A i a c 1 
A1651 


U - z8 -0.25 


1 qc + 1.61 


c OO+0.45 

6 -22l .4i 




cnK+200 
695 -i08 




(k) 


A1664 




1 ro + 1.29 
i - 05 -2.32 












A2029 


6.86+2-S 

— U.oo 


7.03+g-jg 

— U.oz 


7 no+o.39 

— U.OD 


(b) 


1164+ 9 * 

— ( 


93 


(k) 


A2104 




6.85tJ;£ 


Q 1 O + 0.69 


(c) 


1201±200 


51 


(1) 


A2163 


6.63 + i^ 


7 ^ +1 - 78 

' - oo -1.60 


12.3 tH 


(e) 


1698 




(m) 


A2204 


°-' z -1.29 


6.8li^ 


6.38±0.23 


(b) 








A3695 


c c;n+1.47 
°- ou -1.30 


1 7+1-38 
01 ' -1.24 


6.67l?;|| 


(b) 


779t% 


96 


(k) 


A3739 


9 20+ 1 79 


q 41+1.51 












A2345 


5.25_^ 47 


5.93ll;2 












A2384 


^•^-l.OS 


4.02t ; 97 












A2426 










846±100 


15 


(n) 


A3856 


6-37±l;it 








729±142 


22 


(h) 


A3888 


fi 10+1-66 
D.1Z_ 1 4g 


6.46±i;s 


8.46l 3 2 j3 


(b) 


1102±ig? 


50 


(g) 


A3984 


7 kq+1.11 
'■ oy -1.03 


8.30t ;^ 












A2597 


o oq + 1.06 

°- oo -0.93 


4 56 +a82 


4.40±g;| 


(f) 








A4010 


o i 7 +1.00 
01 1 -0.87 


o 00+O.8I 
d.88_ 74 






625l 9 f 


36 


(k) 



Note. — (1) Cluster name. (2), (3) X-ray temperatures inferred from the SIS/SIE 
modeling of the shear data. (4) X-ray temperature with 66% uncertainties from the 
literature. (5) References for Tx- (6) Velocity dispersions found in the literature. 
(7) Number of galaxies used to determine a v . (8) References for <r v . 

References.— (a) Allen (2000) (b) Ikebe et al. (2002) (c) White (2000) (d) 
Valtchanov et al. (2002) (e) Markevitch & Vikhlinin (2001) (f) Markevitch et al. 
(1998) (g) Girardi & Mezzeti (2001) (h) Alonso et al. (1999) (i) Durrct et al. (1998) 
(j) Collins et al. (1995) (k) Fadda et al. (1996) (1) Liang et al. (2000) (m) Squires et 
al. (1997) (n) Mazure et al. (1996) 
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Table 4 

Weighted means and errors of the ratio between cluster temperatures (X-ray and 
dynamical) and weak-lensing inferred values. 





All 


T x < 8.0 keV 


T x > 8.0 keV 






a v < 1122 km s" 1 


a v > 1122 km s" 1 


(1) 


(2) 


(3) 


(4) 


Tx/Tsis 


1.27 ±0.09 


1.14 ±0.03 


1.55 ±0.19 


Tx/Tsie 


1.17 ±0.08 


1.04 ±0.05 


1.42 ±0.15 


(cr v /a S is) 2 


1.20 ±0.13 


1.01 ±0.14 


1.54 ±0.20 


(o-v/ctsie) 2 


1.11 ±0.12 


0.90 ±0.13 


1.41 ±0.15 




5 10 15 



T 3ISiSIE (keV) 

Fig. 6. — Comparison between the ICM tempera- 
tures inferred by the fitting of isothermal profiles 
to the shear data, Tsis and Tsie, and from X- 
ray measurements, Tx- The squares correspond 
to the spherical model and triangles to the ellipti- 
cal. The error bars of the latter were suppressed 
for clarity. The solid line is defined by Tsis,Sie 
= Tx- The dotted and long dashed lines show 
the best obtained by the SIS and SIE models, re- 
spectively, when the origin is kept constant. The 
short dashed line indicate Tx = 8 keV. Clusters 
with higher temperatures show signals of dynamic 
activity. 



(<t v /<tsis,Sie) 2 for both SIS and SIE mass models, 
as well as their errors of the mean. These quanti- 
ties are proportional to the ratio between the mass 
computed from the X-ray emission or cluster dy- 
namics and the mass inferred from weak-lensing. 
For the clusters in our sample for which this anal- 
ysis is possible, we verify that the temperature 
inferred from the lensing analysis is 17% to 27% 
lower than that determined by X-ray observations; 
the squared velocity dispersions determined from 
galaxy velocities are also larger than that obtained 
from lensing, by a factor between 8% and 20%. 
These factors are smaller for the elliptical models. 
Other comparisons between dynamical and weak- 
lensing mass estimates have been made by Irgens 
et al. (2002) and Smail et al. (1997). Whereas the 
latter work has found values of asis ~50% greater 
than a v , the former has found a broad consistency 
between these two quantities. 

Dividing the sample in clusters colder and hot- 
ter than Tx = 8.0 keV (or, equivalently, a v larger 
or smaller than 1122 km s _1 ), the first set presents 
ratios consistent with unity within the errors, 
whereas the second set shows excess (relative to 
lensing results) of 40-50%; see Table 4. Not sur- 
prisingly, all clusters that we have already iden- 
tified as dynamically active belong to this second 
set. 

These results find a natural explanation in 
the hierarchical scenario, where the most massive 
structures are being formed at this very moment. 
If this is the case, higher ICM temperatures are 
expected (compared to an equilibrium state), be- 
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Fig. 7. — Comparison between the velocity dis- 
persion found by fitting isothermal profiles to 
the shear data, ctsie,sie, and those estimated 
through spectroscopic measurements, cr v . The 
squares correspond to the spherical model and tri- 
angles to the elliptical. The error bars of the latter 
were suppressed for clarity. The solid line is de- 
fined by <Jsis,sie = °v are equal. The dotted and 
long dashed lines show the best fit obtained by the 
SIS and SIE models, respectively, when the origin 
is kept constant. The short dashed line indicate 
(T„ = 1122 km s _1 , which correspond the gas tem- 
peratures of Tx = 8 keV. Clusters with higher 
velocity dispersions show signals of dynamic ac- 
tivity. 

cause the gas is heated by shocks produced by 
galaxy groups falling in to the cluster as well as 
mergers with other clusters (e.g. Markevitch et al. 
2002). 

Higher velocity dispersions are also expected, 
since the cluster dynamics will be affected by mo- 
tion streams and substructures, that broaden the 
velocity distribution. 

It is worth investigating further whether Tx ~ 
8 keV may indeed be considered a threshold to 
discriminate between relaxed and active clusters, 
because many studies of the cluster mass function 
are based on X-ray data, that usually select the 
most massive clusters. For example, 11 out of the 



106 cluster of the sample selected by Reiprich & 
Bohringer (2002) have temperatures equal to or 
higher than 8 keV. It is also interesting to note 
that if we use an empirical a — T relation for clus- 
ters as the one of Xue & Wu (2000), we find a 
slightly better agreement between lensing and X- 
ray properties; the main trends discussed before, 
however, remain the same. 

4.3. Comparison of asis with previous es- 
timates 

Some of the clusters studied here have already a 
weak-lensing analysis made by other authors. Ta- 
ble 5 compares our asis with other works that 
used the same model. 

4.3.1. A2744 

The largest difference between the results in Ta- 
ble 5 is for A2744. It has a shear profile very flat 
(see Figure 8). Smail et al. (1997) found for this 
cluster g ~ 0.2 inside a radius of 1', whereas our 
data give g ~ 0.14 in 1.0' < R < 3.2'. Taking 
into account that Smail et al. (1997) have used in 
their analysis galaxies fainter than those consid- 
ered here, thus probably more distant, and also 
closer to the cluster center, it is not unexpected 
the detection of higher gravitational distortions 
by them. Indeed, as discussed in Section 4.4, 
this cluster probably correspond to two structures 
close to the line-of-sight. Isothermal laws are in- 
adequate, and the fact that the xled f° un d for this 
cluster is the highest in this sample confirms it. 

4.3.2. A2104, A2204 and A2345 

There are three clusters studied by Dahle et 
al. (2002) in common with the present sample: 
A2104, A2204 and A2345. For the last two, the 
values of asis obtained by these authors are con- 
sistent with the values presented here, but for 
A2104 there is a large difference. A point that is 
worth mentioning is that these authors have used 
the approximation g ~ 7 instead of g = 7/(1 — k). 
Adopting the same approximation to fit our data 
with the SIS profile, we obtain the following val- 
ues: 1092, 1086 and 957 km s _1 , respectively, in- 
stead of 1039, 1028 and 909 km s -1 . It demon- 
strates that this approximation causes an overes- 
timation of the inferred velocity dispersion by a 
factor of ~5%, depending on the cluster density 
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Table 5 

Comparison of asis with other works 



(1) (2) (3) (4) 

Cluster <Jsis (km s 1 ) crsis (km s 1 ) Rcfs. 
Present work Others work 



A2744 


1491 ± 116 


870l« 


(a) 


A2104 


1039 ± 108 


1390 ± 180 


(b) 


A2163 


1021 ± 146 


740 


(c) 


A2204 


1028 ± 104 


«so:t£8 


(b) 


A2204 


1028 ± 104 


1035+^ 


(d) 


A2345 


909 ± 138 


87011° 


(b) 



Note. — (1) Cluster name. (2) Value of asis 
from the literature. (3) difference between the value 
of asis determined in this work and the value found 
in the literature. (4) Reference of the other weak- 
lensing measurement. 

References. — (a) Smail et al. (1997) (b) Dahlc 
ct al. (2002) (c) Squires et al. (1997) (d) Clowe & 
Shneider (2002) 
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and of the field size, but this is not the cause of 
this difference. In fact, Dahlc et al. (2002) re- 
port problems with their images of A2104. For 
their weak-lensing analysis they used only one I 
band image, obtained with 1 hour integration with 
a seeing (FWHM) of 0.7" at the University of 
Hawaii 2.24 m telescope. Our data for A2104, on 
the other hand, has been collected in one of the 
best nights at VLT in that semester, when the 
seeing was 0.41". This makes us very confident of 
our results. 

Our A2104 mass map is also very different from 
the one obtained by these authors. Accordingly 
to them, the mass distribution is nearly perpen- 
dicular to the major-axis of the central D galaxy, 
which is also along the direction of the two bright- 
est galaxies, and is the direction where the gravi- 
tational arcs are seen. However, our mass map as 
well as the fitting of SIE models show that A2104 
has an elliptical mass distribution aligned with the 
central galaxy major-axis, as expected if the BCG 
light is a good tracer of the cluster potential (see 
for example Mellier, Fort & Kneib 1993). 

4.3.3. A2163 

Our estimate of the velocity dispersion of A2163 
is also much higher than that reported by Squires 
et al. (1997). These authors have observed this 
cluster in V and I with total expositions of one 
hour at the 3.6 CFH telescope, under seeing of 
~0.8". They have used ~ 700 galaxies fainter than 
1=20.5 mag or V=22.0 mag for the lensing analy- 
sis, in a field of 7 x 7 arcmin. 

Here, we have used a much smaller sample of 
background galaxies, in a field of comparable size, 
because we have considered only galaxies fainter 
than 23.3 R mag, that is roughly 2.0 mag fainter 
than the bright end limit used by Squires and col- 
laborators. The images of this cluster were taken 
in the same night as those of A2104, when the see- 
ing was exceptional. A2163 is at a moderately high 
rcdshift (z = 0.2) and lies in a region of the sky not 
too far from the galactic plane (b=30.5°), where 
the star field is dense and the galactic absorption is 
not negligible (Ay = 1.74; Schlegel, Finkbeiner & 
Davis 1998), so the weak-lensing sample of Squires 
et al. (1997) may be more affected by contamina- 
tion by stars and foreground and cluster galaxies 
than the present one. This might be the cause for a 
possible dilution of the lensing signal measured by 



them. Indeed, adopting a comparable magnitude 
cutoff, R > 21.5 mag, we obtain asis =887±145 
km s _1 , that is much closer to the Squires et al. 
(1997) result. 

4.3.4. A2204 

Table 5 shows that our results for A2204 are 
in very close agreement with those from Clowe & 
Shncider (2002). This work is based on images 
taken with the wide-field camera of the ESO /MPI 
2.2m telescope (34' x 34'), indicating that when 
the data are well described by an isothermal mass 
profile, as is the case of A2204, the small field 
of the FORS1 camera does not introduce relevant 
systematic errors. 

4.4. The peculiar case of A2744 

This cluster, also known as AC118, is inter- 
esting for several reasons. It is the most distant 
(z = 0.308) and apparently the most massive clus- 
ter in our sample. It is also the only cluster where 
we have found gravitational radial arcs (see Cam- 
pusano et al. 2003), and has a flat shear profile 
which is very poorly fitted by isothermal mod- 
els. Furthermore, its weak-lensing mass is smaller 
than the virial mass but, surprisingly, significantly 
larger than the X-ray estimate. We should have 
in mind that the use of an isothermal profile to 
fit a shallow mass distribution tends to underes- 
timate the cluster mass, since the shear depends 
more on the density gradient than on the density 
actual value. But an important clue for the inter- 
pretation of this cluster comes from the dynamical 
analysis of Girardi & Mczzeti (2001), who found a 
velocity distribution with two superposed peaks, 
produced by two structures along the line of sight, 
with velocity dispersions of 1121 and 682 km s _1 . 
Weak-lensing, in this case, is sensitive to the to- 
tal mass, but single isothermal mass models are 
clearly inadequate for this system. On the other 
side, the mass inferred from the X-ray analysis will 
be dominated by the luminous and hotter struc- 
ture; for a v = 1121 km s _1 a value of Tx = 8.0 
keV is expected, somewhat smaller than the X-ray 
measurement (see Table 3). In this scenario the X- 
ray luminosity is the sum of the luminosities of 
both structures. The value of Lx for A2744 is ac- 
tually very high (22.05 x 10 44 erg s _1 ), ranked in 
second place in this sample. Of course, both Tx 
and Lx may be increased with respect to their 
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equilibrium values by shocks between these two 
substructures and, therefore, all these results are 
consistent with the hypothesis that A2744 is actu- 
ally two close structures along the line of sight or 
even a pair of clusters in process of merging. 

4.5. Negative weak-lensing detections: A1651 
and A1664 

For only 2 out of 24 clusters (8.3%) we were 
unable to reconstruct the two-dimensional mass 
distribution. This number is comparable to that 
obtained by Dahle et al. (2002), who failed to de- 
tect weak-lensing in 1 cluster out of 39 (2.5%). 
These authors have adopted a sample threshold 
of Lx = 10 45 erg s _1 , two times larger than ours. 
Our two failure cases are clusters with luminosities 
below this threshold. 

For A1664 we were unable to reconstruct the 
mass density distribution and to fit a SIS profile, 
and we have got a very marginal detection when 
fitting a SIE profile. This cluster is at z = 0.13 and 
has a X-ray luminosity very close to our threshold. 
Neither the X-ray temperature nor the velocity 
dispersion is available in literature for this object. 
Moreover only 246 galaxies were used for the weak 
lensing analysis (the third smaller number of the 
sample). It seems that we were unable to detect 
the weak-shear because its signal is small and is 
dominated by the noise in its distortion map. 

The case of A1651 is less clear, since this cluster 
seems to be more massive than A1664 (see Tabic 
3). In fact, the SIS fit resulted in a very small value 
for a sis, but the SIE fit led to a significant detec- 
tion, resulting in a value of asis = 539 ± 198 km 
s _1 , consistent with the dynamical velocity disper- 
sion found in literature, 695^Qgkm s_1 (Fadda et 
al. 1996). The position angle of the SIE model is 
almost coincident with that of the BCG major- 
axis, as found for most of the clusters in our 
sample. However, the lensing signal is still too 
low to allow a successful mass reconstruction by 
LensEnt. 

Another common factor for these two cases 
is that seeing conditions during the observations 
were within the worst quartile of the sample. This 
puts in evidence how critical seeing conditions are 
for ground-based, weak lensing studies. 



5. Conclusions 

We analyzed the mass content of galaxy clusters 
belonging to a well defined sample of 24 Abell clus- 
ters brighter than h x =5 x 10 44 erg s" 1 (0.1-2.4 
keV Ebeling et al. 1996), spanning over the red- 
shift range 0.05-0.31, using current techniques of 
weak- lensing analysis and homogeneous observ- 
ing material of subarsccond image quality. The 
resulting catalog of mass maps determined for 22 
of these clusters, together with the correspond- 
ing light maps, have been put together in the Ap- 
pendix. Our main conclusions may be summarized 
as follows: 

1. We were able to detect significant weak- 
lensing signal in 22 out of 24 clusters. This 
high success rate shows the feasibility of 
weak-lensing studies with 8m-class tele- 
scopes using service mode observations and 
relatively short exposure times. It also in- 
dicates that the X-ray luminosity is indeed 
a good way to select massive clusters. Non- 
detections in A1651 and A1664 are probably 
due to a combination of poorer observing 
conditions with low mass content. 

2. The center of the mass and light distribu- 
tions of the clusters are coincident for ~77% 
of the sample (17 out of 22). 

3. Few clusters present massive substructures, 
what can be due, in part, to small fields (0.4 
to 1.8 Mpc on a side) and the relatively low 
resolution of the mass maps. When signifi- 
cant substructures are seen, they are gener- 
ally associated with bright cluster members. 
However, at least in one cluster (A1451) 
there seems to exist a substructure without 
a clear optical counterpart. 

4. The clusters analyzed here present impor- 
tant departures from spherical symmetry, as 
can be verified by the better fits obtained 
with elliptical profiles. 

5. We have found, for the first time in an sta- 
tistically significant sample, that the dark 
matter and brightest cluster galaxy major- 
axes are strongly aligned: for 62.5% of the 
clusters (15 out of 24) the difference between 
their position angle is smaller than 20°. 
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6. Most clusters are in or near a state of dy- making their softwares. im2shape and LensEnt 
namical equilibrium. This diagnosis derives available to us. 

from the agreement between their velocity 
dispersions and the temperature of their 
ICM, directly measured and/or inferred 
from weak-lensing data. Except for A2744, 
A1451 and A2163, which also present evi- 
dence of substructures or other complexities, 
the other clusters show agreement between 
these quantities at a 1.5 a level. 

7. Clusters in our sample with Tx > 8 keV (or 
a v > 1120 km s _1 ) show signal of dynami- 
cal activity. A2163 and A1451 present large 
differences between lensing and dynamical 
mass estimates and seem to be far from equi- 
librium. In both cases, this conclusion cor- 
roborates previous X-ray analysis. 

8. Abell 2744 is the single cluster in this sample 
that has X-ray measured temperature more 
than 1 a below the value inferred by weak- 
lensing. Taking also into account the com- 
plex dynamics of this cluster, we may explain 
this discrepancy by assuming that A2744 is a 
superposition of two clusters along the line of 
sight, near each other or in process of merg- 
ing. 

Most of these conclusions supports a hierarchi- 
cal scenario where massive bodies are formed by 
the agglomeration of smaller ones, and the depar- 
tures of equilibrium described above are indeed an 
evidence that some clusters that are at the top of 
the mass function are still in the process of active 
evolution. 

This work is based on observations collected at 
the European Southern Observatory, Chile (ESO 
P67.A-0597). ESC and LS acknowledges sup- 
port by by Brazilian agencies FAPESP, CNPq and 
PRONEX, and ESC is grateful to the Latin Amer- 
ican Astronomy Network/UNESCO for funding 
traveling to Chile. JPK acknowledges support 
from CNRS and Caltcch. ESC also thanks the 
hospitality of the Department of Astronomy of 
the the University of Chile and of the Labo- 
ratoire d'Astrophysique of the the Observatoire 
Midi-Pyrenees. We are also grateful to Gastao 
Lima Neto and Ronaldo de Souza for useful dis- 
cussions and Sarah Bridle and Phil Marshall for 
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A. Mass and Light Maps 

In this appendix we present figures referent to the weak-lensing analysis. In each figure there are 4 panels. 
In the first panel (from left to right) we show the cluster image with superimposed sticks that represent the 
average ellipticitics and position angles of the background galaxies. The size of the stick is proportional to 
the ellipticity. The horizontal stick in the upper right corner of this panel represents an ellipticity of 0.1. 
We also present an ellipse with the axial ratio and position angle of the best fit SIE profile. In the second 
panel, we show the mass maps, that is, the surface density distribution reconstructed from the shear data, in 
unities of the critical surface density. The contour are in intervals of 8 x 10 11 M Q Mpc~ 2 . In the third panel 
the distribution of light coming from cluster members that belong to the cluster red sequence, in arbitrary 
unities, smoothed with a Gaussian similar to the one used to produce the mass map. The countour are in 
intervals of 1 x 10 11 L Q Mpc -2 in the V band. The fourth panel presents the radial shear profile, where the 
filled circles represent the average cllipticities projected tangentially to the direction to the cluster center, 
and the open squares show the cllipticities project in a direction 45° of the tangential one. The solid line 
represents the best-fit SIS profile. 



A version of this paper with all figures can be found at 
http: / /www. astro. iag.usp.br/~eduardo/shear.ps.gz 
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